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Detailed protocol 


Protocol for co-expression analysis in Figure S1 of Research Article, 

“Controlled hydroxylations of diterpenoids allow for plant chemical defense without autotoxicity” 

DOK: 10.1126/science.abe4713 

Figure S1A. 

The correlation between known genes was developed using package “ggplot2”, with the script attached as “Fig. S1A_heatmap_ggplot2.R” 
Figure S1B. 

The heatmap in B panel was developed using the original script of “heamap.2” in package “gplots”, attached as “Fig. 
S1B_heatmap_gplots.R’. 

Figure S1C. 

The dendrogram was produced using the modified R package “WWGCNA’”. The script was attached with name as “Fig. S1C_Module detection _WGCNA.R’. 
The heatmap in Figure S1C was produced using the exact same script with “Fig. S1B_heatmap_gplots.R’. 


note that since .R files are not an allowable format for this site, all three .R files have been converted to .txt files 
Related files 

Fig. SLA_heatmap_ggplot2.txt 

Fig. S1B_heatmap_gplots.txt 


Fig. S1C_Module detection. WGCNA.txt 
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